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INTRODUCTION 

The subject of geometric modeling is receiving considerable attention 
from the industrial community not only within the United States, 
but throughout the v/orld. Current thinking points toward the 
establishment of an application-independent geometric model within 
a data base framework. The modr;l can then be accessed and used 
by any number of technologies within both engineering and manufacturing. 
This basic model may be used to create application-specific models 
such as a Nastran finite element model or a model appropriately 
paneled for aerodynamic analysis. Such a model may also be used 
by manufacturing to support their numerically controlled operations. 

Thus far a universally acceptable geometric model has not been 
found. However, certain attributes appear to be getting strong 
backing. First, the model should support the use of Boolean 
operators, i.e., UNION, INTERSECTION, and DIFFERENCE. A second 
property of the model is that the bounding surfaces be described 
parametrically. The detailed m.athematical form of the surface 
representation has yet to be resolved. 

The use of the Boolean operator's is attractive since it allov/s the 
user to construct a complex mooel by appropriately combining a 
series of simple models. The use of these operators leads to the 
concept of implicitly and explicitly defined surfaces. With an 
explicitly defined model, the surface area may be computed by 
simply summing the surface areas of the bounding surfaces. For 
an implicitly defined model, the surface area computation must 
deal with active and inactive regions. To illustrate this, 
consider a simple wing-body combination.. To model the combination 
explicitly, the modeler must first capture the wing-body inter- 
section line which will act as a boundary curve for both the wing 
and the fuselage. Since the two components are generally designed 
separately, it would be advantageous to model them separately. 

Requiring the intersection curve to serve as the surface boundary 
curve would probably lead to local distortion in both the wing 
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and the fuselage. In the implicit modeling scheme, both components 
are modeled separately with no consideration of the intersection 
curve. When combined in their desired orientation, the intersection 
curve can be determined numerically. The intersection curve serves 
as the boundary betv/een active and inactive regions. For the wing, 
the portion contained within the fuselage is considered inactive. 

For the fuselage the hole made by the wing is inactive. Thus, 
to compute the surface area of the combination, the area of the 
inactive regions must be disregarded. The advantage of the 
implicit system can be appreciated if one considers the simple 
exercise of moving the wing forward by an inch after the mode! 
has been created. With an explicit model, the entire job must 
be redone, while in the case of an implicit model, only the 
intersection curve need be regenerated (by the computer, not 
the user). 

Examination of the seftv;are required to support implicitly defined 
models identifies the surface intersection capability as a critical 
item. The purpose of this report is to describe a particular 
solution to the surface intersection problem. Included will be an 
overall description of the strategy employed plus a detailed 
presentation of the analytical and numerical considerations required 
to implement it. Finally, a description of the resulting software 
will be given, including user information and examples. While 
the solution is implemented for the specific case of parametric 
cubic surface representation, it is felt that the solution 
should be equally applicable to other mathematical forms. 


2.G SURFACE INTERSECTIONS - AM OVERVIEW 

The purpose of this section is to describe a computational .scheme 
for obtaining the intersection between tv;o surfaces. Sinca this 
problem is quite difficult and the required analysis quite complicated, 
it was decided to begin the discussion with an overview of the strategy 
as well as the definition of terms to be used later in the detailed 
breakdown of the various computational phases. To begin with, the 
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term surface needs to be more accurately defined. VJithin the con- 
text of the present problem, a surface S is defined to be a rectangular 
mesh of parametrically defined patches which approximate or represent 
a surface in physical space. Figure 2.1 shovrs a typical surface 
composed of a 4 x 5 mesh of patches. The figure also indicates a patch 
numbering convention v/hich is important since the patches must be 
stored and retrieved in some logical order. The numbering convention 
also establishes connectivity between the individual patches. An 
additional convention which is not obvious from the figure is that 
the vector product betv/een the u-derivative and the v-derivative 
results in an outward surface normal. 


_ _ _ 


17 

18 

19 

20 


13 

14 

15 

16 

9 

10 

11 

.12 

5 

6 

7 

8 

1 

2 

3 

4 

0 

2 

\ I 

4 ' 


FlflURE 2.1 - A Typical Surface 


Nu = 4 
Nv = 5 


Consider now the problem of finding the intersection between two 
surfaces, Si(u,v) and $ 2 ( 5 , t). The calculation may be decomposed 
into three distinct phases: the h unting phase , the tracing, phase, , 

and the sorting phase . The overall strategy is necessitated by 
the requirement that the calculation be capable of handling 
relatively complex surfaces. The hunting phase locates discrete 
starting points required for the curve tracing operation. The 
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tracing phase creates strings of points lying on the surface inter- 
section. Finally, the sorting phase orders the point strings and 
separates them into disjoint segments or loops. 

2.1 Hunting Phase 

The purpose of the hunting phase is to locate a number of points 
on the intersection v/hich serve as starting points for the 
tracing phase. At least one point on each loop of the inter- 
section curve must be found during this phase, or that particular 
loop will go undetected. This type of situation suggests the 
need for some sort of refinement mechanism which can be made 
to adjust to the most complex of cases. The capacity for 
refinement is obtained by introducing the hunting grid on the 
primary surface Sj. Figure 2.2 shows a typical primary surface 
with such a hunting grid superimposed. By specifying a single 
integer Nq, a grid is defined on Si. The example intersection 
is composed of two loops. It would appear that a smaller value 
of Nq could have been used without losing either of the loops. 
However, if a sufficiently small Nq (Nq = 3} were chosen, the 
loop located in the upper left corner could escape detection. 

The rectangular regions defined by the hunting grid are numbered 
in the same manner as are the conventional surface patches 
(indicated by circled numbers in Figure 2.2). The hunting 
phase consists of taking each curve defined by the grid and 
obtaining its intersection with S^. By taking first the curves - 
of constant u and then the curves of const <it v the intersection 
points are located in a sequence indicated by their numbering 
in Figure 2.2. 

The problem of determini no' the intersection of a given curve 
v/ith $2 requires additional comnent. As in the overall problem, 
the major concern here is to not miss any intei"section points. 
This problem is somewhat anc'.logous to determining the roots 
of a polynomial. The scheme must be able to isolate the 
roots and then determine them. The approach used to isolate 

MGOOMiWELL OOUOt.Aa 

4 



MDC J7789 


the intersection points is a refinement of the hunting grid 
philosophy. The curve is segmented by selecting an integer 
Mq whose role is analogous to that of Ng in the hunting grid. 
Mg should be sufficiently large that two intersection points 
do not occur v/ithin the same segment. In keeping with the 
polynomial analog, a function is sought that will identify the 
presence of an intersection point by exhibiting a sign change 
between consecutive nodes. Consider the vector function: 

V = Pc - Ps = (F.G.H) 

which represents the directed distance between a point on S 2 
and a point on the curve C. For an intersection to occur 
betv^een nodes i and i + 1 on C, it is necessary that 



FIGURE 2,2 - Example of Hunting Grid (Ng = 6) 
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The marching procedure is illustrated in Figure 2.3 



FIGURE 2.3 - Marching Procedure for Curve/Surface Intersection 

This seemingly simple concept requires one final assumption. 

A point on the curve ?c is compared to the point on the surface 
Ps which is closest to ^c. Thus, the vector f has three 
components in parametric space: the curve parameter ( u or 

v) and the two surface parameters s and t. Since the 
effectiveness of this algorithm is not dependent on the 
precise determination of s and t such that |vj in a minimum, 
an approximate scheme may be used. Specifying a value of the 
curve parameter (u or v) determines a point ?*. The use of 
th^ surface inversion algorithm for a small number of iterations 
will satisfy the requirement for the approximate closest 
distance from a point to the surface ?s(s ,tl-?’*=Q. Having 
isolated an intersection point, 

i.e. (u or v)^ £ (u or v) _< (u or 

the full three- variable problem (’^=0) is then solvec subject 
to the indicated constraint on the curve parameter. The hunting 
phase will produce a table of intersection points lying on the 
grid boundaries of Si. The ordering of this table obviously 
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depends on the order in which the grid curves are taken as 
well as the direction of march along these curves. To 
illustrate this, the grid intersection points are numbered 
in Figure 2.2 as they would appe?r in the table. 

2.2 Tracing Phase 

The hunting phase produced a number of discrete inter- 
section points lying on the boundaries of the hunting 
grid. As shown in the example (Figure 2.2) the ordering 
of these points is not optimum regarding the further 
enrichment of the data to form the complete intersection. 
The intersection curves are traced across the interior 
regions of the hunting grid in the following manner. 

First the regions are analyzed one at a time. The 
ordering is logically similar to the patch numbering 
scheme and is indicated in Figure 2.2 by circling the 
region numbers. In eacv, region, the boundaries are 
compared with the tabulated points. When a match is 
found, the portion of the particular loop containing 
that boundary point and lying within that region is 
traced, usually ending at another boundary point. Thus, 
the direction of trace across a region is established 
by which endpoint occurs first in the table. The ending 
point is excluded from the rest of the search ( for that 
region ) in order to prevent multiple tracing of the 
same segment. In this manner, regions containing more 
than one intersection segment (more than tv/o boundary 
points) will be correctly handled. 

At the starting point (for a given region) the tangent 
vector to the intersection curve may be obtained by 
taking the vector product betv/een the normals to surfaces 
Si and S 2 . The order of the vector product is established 
such that the tangent vector points toward the interior of 
the region. This tangent vector provides the deriva- 
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tives of the four parameters (u, v; s, t) with respect to arc 
length along the intersection curve. Using a step siz* based 
on curvature and its relation to slope change along the inter- 
section curve, all four parameters are incremented, e.g. 

“in “ “1 * (^), '‘“1 

This step places the next point in the neighborhood of the 
intersection curve and a secondary iterative convergence is 
required to achieve the desired tolerance. Since there are 
four unknowns (u, v, s, t) and only three equations, one of 
the unknowns must be held constant. By selecting one of the 
four possible groupings of unknowns, the system is made 
determinate and a converged solution obtained quite rapidly 
using Newton's method. This two level marching process is 
continued until a region boundary is encountered. 

2.3 Ordering Phase 

The example shown in Figure 2.2 illustrates the need to sort 
and order the various point strings into loops. The order 
in which the segment strings v/ere produced is indicated by 
the lower case letters. While this operation might be 
regarded as mere bookkeeping, it is significant enough to 
warrant discussion. The result of performing the first two 
phases is a collection of N point strings (N = 10 in 
Figure 2.2). With regard do sorting; the only significant 
points are the first and last points in each segment string. 
Beginning with the first string, the start and end point of 
the current loop (containing that string) are compared to 
the remaining strings. When a match is found, that string 
number is saved in an ordering array as a signed integer. 

The minus sign is used to indicate that a particular string 
requires reversal when it is concatenated to form a loop. 

In the example shown in Figure 2.2, the ordering arrays are 

0€SiJGt.A9 ATIGM 
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as follows: 

LOOP 1 : a, -b,-c, d, -e,-f 

LOOP 2 : -h, g, -1", j 

Some further explanation is required regarding LOOP 2. When 
LOOP 1 is completed, only four segments remain: g, h, i, and 
k in that order. The scan to order LOOP 2 begins with g 
which has boundary point 1 as its beginning and point 10 as its 
end. The scan proceeds to segment h where a match is found 
with the start point of LOOP 2. Since the direction of LOOP 2 
has already been established by g, h must be reversed and 
should precede g in the ordering table. At this stage, 

LOOP 2 has point H as its start and point 11 as its end. 

The top end of i matches the top end of LOOP 2 so it must be 
reversed and point 2 becomes the top end of LOOP 2. Finally 
the bottom of k matches the top of LOOP 2. At this stage, 
there are no segments not accounted for and the sorting phase 
is complete. 

One final consideration of the sorting phase is illustrated in 
Figure 2.4. LOOP 1 is started by a running from points 2 to 
4. When b is compared to LOOP 1 it is found to be disjoint. 

At the end of the first pass, the ordering array for LOOP 1 
is a, -c, d, v/ith b considered as disjoint. This failure is 
a result of the order of the search and is easily remedied by 
performing a second pass prior to starting on LOOP 2 and would 
append -b to the ordering array for LOOP 1. 



FIGURE 2.4 - Sorting Example 
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3.0 SURFACE INTERSECTIONS - MATHEMATICAL DETAILS 

Most of the mathematical manipulations involved with solving the 
surface intersection problem concern the solution of three simultaneous 
equations - one for each of the Cartesian coordinates of our three- 
dimension physical space. The distinguishing feature in these 
subproblems is usually the number of unknown variables. The general 
surface intersection problem has four unknowns - the parametric 
coordinates of each surface {u, v: s, t}. Since the solution is a 
space curve, the solution actually involves determining these four 
unknowns as a function of single parameter, e.g. the arc length o. 

This problem might be contrasted to that of determining the point 
on a space curve closest to a given point in space. The latter 
problem involves only one unknov/n - the parametric coordinate of the 
point on the curve. The following subsections will identify the 
various subproblems occurring in the general intersection problem 
and provide the mathematical basis for their solution. 

3.1 Determination of Parametric Derivatives from a Tangent Vector 

It was stated in Section l.Z that the vector product of 
the two surface normals at a point on the intersection 
curve produces a vector tangent to that curve. Furthermore, 
this tangent vector may be used to obtain the derivatives 
of the parametric coordinates with respect to arc length 
along the intersection curve. Let the tangent vector t 
be given by 




T = Ni X No = — lx— 1 X — 2x^i2 
' 1^1 X 1^2 3U X 3s ^3t 

and T s 7/\l\. 

A 

The unit vector T satisfies the equation 


A dP 
T = do 


MCDOMMELL, OOUGt-AS COKfOEATIOM 
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which for $i may be expanded as 




where the primes are^ used to denote derivatives with respect 
to arc length a. A similar expression can be written for S 2 
in terms of s^ and t^. 

Each of these vector equations represents a set of three 
nonlinear algebraic equations in two unknowns. Thus, one 
would expect a degree of freedom in their solution, i.e., 
any set of two equations should suffice. By taking the 
appropriate vector piroducts, u^ and v^ may be isolated as 
follows 

U‘ » - X f/ili 


where Ni E X |E 

v> = X T/N, 

Some explanation is required for the above equation in order 
to properly interpret the quotient of two vectors as a scalar. 
The equations are evaluated by selecting any of the three 
vector components which is where the added degree of freedom 
comes in; e.g. 


/8y dz\ 

\3V 3o/ 

-1 

lav dg/ 

/3y 

\3U 3vy 

> /ay 3z\ 
' lav au/ 


where the x-component was selected. Barring the presence of 
more important considerations, the component of W having 
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the largest magnitude will be selected. 

3.2 Surface Inversion 

The problem of inverting the surface eouations to solve for 
the parametric coordinates of a given point in physical 
space is similar to the previously discussed parametric 
derivative determination in that it is also overspecified. 
The set of vector equations v/hich must be solved is given 
by 


V = P(u,v)-P* 


0 . 


three equations in tivo unknowns. Since these equations 
may not be inverted analytically, an iterative scheme is 
required. The first variation of the above equation may 
be written as 

* I? 


Given an initial guess (ui, vi ) such that V(ui,vi)=Vi, 
the variations in the above equation may be interpretted 
as 

n+i - V, ={|£)i (ui+1 - Ui) (n^., - Vi) 


Setting Vi+i = 0 and taking the appropriate vector products, 


Ui+i = Ui + Ai/Ni 


Vi 4-1 = Vi - Bi/Ni 

Ai= V, 


H- fix Vi 
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In the above expressions, the ratio of tv-;o vectors is inter- 
pretted in the same way as it was in Section 3.1. If P* lies 
on the surface, the solution will tend to converge to the 
point on the surface closest to In the latter case, the 
cutoff criterion may be as simple as exceeding a fixed 
number of iterations or possibly 

IVi + i - Vi I <€ 

3.3 Curve-Surface Intersection 

The cuwe-surface intersection problem is the basis for the 
hunting p’lase. It differs from the above two problems in 
that it results in a fully specified system of three equations 
in three unknowns. Let the curve be given as Pc (v) and 
the surface by P<; (s,t). The system of equations is given 
by 


V = Pc (v) - Ps (s,t) = 0 
The first variation of this equation is written as 


6V 





5t 


The variations in the parameters may be isolated by taking 
appropriate vector and scalar products. First, for example 












Furthermore 


i£c 

3V 



X 6V 
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In a like manner, the iteration is supported by the 
equations 


= Vi - 

fl£si 

1 3S * 

\ at 

xvO' 

l/D 

^i+1 ' 

1 3Psi 
1 at 

/3Pci 
\ av 

X Vj)] 

/D 

ti = ti + 

[i!ci 
1 av 

/iEsi 
\ as 

X Vl) 

I/D 


where 





It should be noted that since this system is fully specified, 
the redundancy as characterized by the ratio of vectors is 
no longer present. 

Surface- Surface Intersection 

In Section 2.2 it was implied that final convergence in the 
tracing phase involved an underspecified system, namely three 
equations in four unknowns. 

^ = Pi(u,v) - Pz(s,t) = 0 


It was also stated that this problem v;as resolved by holding 
one of the unknowns constant. This would reduce the problem 
to a curve-surface intersection already covered in Section 
3.3. The criterion used to determine which variable to fix 
is of some interest. The most obvious criterion would be 
to fix the variable which had the smallest derivative with 
respect to arc length at the previous converged point. A 
second approach begins by observing the similarity between 
the overspecified and the underspecified problem. In both 
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cases a decision must be made before the computation may 
proceed. In the overspecified problem the decision was 
based on maximizing the magnitude of the denominator. By 
observing that any decision in the underspecified problem 
reduces 1t to a fully specified problem whose solution also 
Involves a denominator, the same criterion could be applied. 
The variable is fixed that provides the maximum denominator 
in the resulting fully specified problem. This criterion is 
primarily an intuitive one with the only obvious advantage 
being avoidance of numerical problems associated with small 
1 d[. The resulting four sets of equations are as follows. 

I. u a constant: Uj^^ = u. 
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IV. t = constant : = Ui + 


3fii • |P2i^ ^ )j 

av \ as ^ 


'i+1 




^i+i = Si - 


1^- (|^x^i)l/Dt 


t 


i+1 


Dt = 


= ti 


->• -¥ -> 



3.5 Step Size Seleetion for Tracing Phase 

At a point on the intersection curve, the curvature may 
be expressed in terms of the angular rate of change. 
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dc 



FIGURE 3,1 - Step Size Selection 


By selecting a nominal value of the angular change, the 

step size 6a can be obtained from 6a = 4^. The tangent 

^ ^ ^ ^ 

vector to the intesection curve T is defined by T = eNj x N 2 » 
where 

Ni = 1^ X = a vector normal to Si, 

No = 

^ ^ ^ = a vector normal to Sj, 


e = + 1 = a scalar multiplier to correct the sense of 
the tangent vector. 

The curvature jK| may be defined as |-^j where t = f/|tj. 
Expanding this expression, 


df ^ 4_ f dt 

da |T| I da 


(T 




Since T is a function of all four independent variables 
(u, V, s, t). 



The parametric derivatives (u^, v^, s^, t^) are derivable from 
f (see Section 3.1). Expanding this equation yields 

MGOOMtVELL. OOUCH-AE ATIOM 
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Thus, the step size may be estimated based entirely on t<ie 
local surface properties. 

.6 Interpolation of the Intersection Point String 

The surface intersection computation yields an ordered 
string of points lying on the intersection curve Cwithin a 
given tolerance). The string is subdivided into closed 
curves or loops. The problem to be considered here is that 
of interpolating this output data and it is sufficient for 
this purpose to consider a single loop. The number of 
output points is controlled by the spacing constraints on 
maximum step size and maximum slope change. The values of 
these constraining parcuiieters depend heavily on the use 
intended for the output data. Although a comprehensive 
study regarding the use of the output data in a variety 
of applications is beyond the scope of the present effort, 
it is felt that the issue should be discussed. 


Let us begin with the premise that most computer programs 
are somewhat limited by available storage. The intersection 
calculation is capable of producing a large amount of output 

ntcoo/yiyeut. ooifGLAa eo#7#>o«v/irfOM 

18 



MOC J7789 


data. Consider the case of a single loop intersection that 
produces 100 output points. A single point consists of the 
following data: 

1) (x,y,z) - physical coordinates of the point 

2) (u,v,s,t) - parametric coordinates of the point 

■ derivatives of the parametric 
coordinates with respect to arc length along the 
intersection curve. 

Thus the example case provides 1100 words of data. Obviously 
some of the data is redundant and could be derived from the 
other data. There would appear to be a good argument for 
tailoring the output to fit the application. If a standard, 
nonredundant output were used, e.g. the physical coordinates, 
applications which required the other data would have to pay 
the computational penalty. 

Regardless of the particular set of data selected, a decision 
must be made concerning its output format. In particular, 
should the data be returned as discrete points or should it 
be curve fit? Once again, the application will probably 
provide the answer. Consider the example case where it has 
been decided to work with the intersection in physical space. 
The discrete point format would provide 300 words of data. 

An PC fit of the data, using 10 segments in order to satisfy 
some tolerance constraint, would produce 120 words of data. 
Clearly some compromise must be made betv/een core requirements, 
execution speed and accuracy. 

If the application calls for the parametric coordinates there 
are additional options available if curve fits are desired. 

All four parametric variable could be defined as functions of 
a fifth parametric variable monotonically related to arc 
length along the intersection curve. In this case, the 

MCOOMMELL DOUGUAS COG^GEATiOM 
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tolerandng can be carried out in physical space with the 
error at any given point given as 

6 = iP(u.v) - Q(s,t)| 

where (u,v,s,t) are taken at the same value of the independent 
parameter. A second approach was to treat each of the pairs, (u,v) 
and (s,t), as separate problems to be fit independently. The entire 
fitting and tolerancing is carried out in the parametric plane so 
that the tolerancing loses its relationship with the physical inter- 
section curve. The first approach provides interpolation with 
physically meaningful precision, but generally requires a large 
number of segments resulting in a large quantity of data. The 
second approach sacrifices the controlled precision but significantly 
reduces the amount of data which must be saved. The question as 
to which approach is best is highly anplicatlon dependent and really 
needs further investioation. The second approach is currently 
incorporated in the subject program. 

SOFTWARE DESCRIPTION ‘ 

The purpose of this section is to describe the computer program 
written to implement the strategies outlined in the previous sections. 

The code was written in FORTRAN IV for implementation on the IBM 
370 timesharing system, TSO. The input to the program has two forms. 

First the actual coefficient data for the tvjo surfaces to be operated 
upon is supplied in the form of formatted input files which are 
automatically read in by the main program. The second group of data 
constitutes a set of prompted input parameters to be entered from the 
terminal by the user. The second Set of data consists mostly of 
tolerances, step sizes, or option flags which are likely to vary 
from case to case. The program output is written to an output 
file whose later disposition is a user option. 

4.1 TSO Command List 

The use of a CLIST attempts to remove the IBM job control 
language burden from the user by submerging it. The purpose 

MeoOMMEt-L. OOC/<Sfb^9 COEfOnATBOM 
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of the CLIST in the present case is to handle the allocation 
of the various files used by the program during execution. 

The CLIST has two positional parameters plus an optional 
keyword parameter. The CLIST is invoked by the following 
statement: 

EXECUTE NAME 'INI INL OUTPUT(NAMEOUT) ' 

where NAME is the name given to the CLIST, INI and IN2 are 
the names of the data sets containing the coefficient data 
for the ti/io surfaces, and NAMEOUT is the name to be given 
to the output data set. If the optioi-al keyword parameter 
is not specified, the output data set is given the default 
name OUT. DATA. For completeness, the CLIST is as follows. 

PROC 2 INI IN2 OUTPUT(OUT) 

FREEALL 

ATTR FI BUFNO(l) 

ALLOC F(FTOIFOOI) DA(&IN1 , .DATA) SHR USINR(Fl) 

ATTR F2 BUFNO(l) 

ALLOC F(FT02F001) DA{&IN2. .DATA) SHR USING(F2) 

ATTR F3 BUFNO(l) 

ALLOC F(FT03F001) USING(F3) NEW BL(3200) SP(20,10) 

ATTR F4 BUFNO(l) REC{F B) LR(140) BL{1680) 

ALLOC F(FT04F0Q1) DA(&OUTPUT..DATA) USINR{F4) MEW BL(1680) SP(10,10) 
ALLOC F(FT05F001) DA(*) 

ALLOC F{FT06F001) DA{*) 

CALL NAME LOAD(TEMPNAME) 

FREE F(FTOIFOOI) 

FREE F(FT02F001) 

FREE F(FT03F001) 

FREE F(FT04F001) 

END 
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4.2 Surface Data Files 

The surface data consists primarily of parametric bicubic 
patch coefficients in geometric form. The patches are 
ordered as indicated in Figure 2.1. For completeness, a 
parametric bicubic patch may be represented in geometric 
form by a matrix of vectors B which may be written as 


p(o.o) %j,i) |£(o.o) l^(o.i) 
p(i.o) p(M) I^O.o) l^d.D 
f„(0.0) |f(0.1) 3|-„(0.0) 3^C0.1) 




a2p 

auav 


( 1 , 0 ) 


a ^ 
auav 


0 . 1 ) 


The first line of data in the surface data file contains 
the following information 


NP = number of patches 

ANAM(3) = three word hollerith title 

NU = number of patch columns (see Figure 2.1) 

NV = number of patch rov'S (see Figure 2.1) 

The second block of information contained in the surface 
data file is a transformation matrix which is used to 
orient the surface in space. The matrix is a 3x4 array 
which may be partitioned into a 3x3 rotation matrix R and 
a 3x1 translation vector T. Since the geometric form of 
the patch coefficients consists of vectors, the surface may 
be transformed by trar«forming the defining vectors. The 

f 

only notev/orthy feature of implementing the transformation 
is that all vectors are rotated, but only the corner point 

position vectors are translated. The transformation matrix 

is v/ritten using FORTRAH (3D20.13). 
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The remainder of the data consist*; of actual patch coefficient 
data. Each card image is v/ritten using FORMAT (3D20.13,2I5,6x, 
A4), Each card image contains one of the vector elements of 
the B matrix. The tv;o integers indicate the element number 
and the patch nimber, vfhile the hollerith text maybe used for 
identification purposes. The element number varies fron 
1 through 16 listing the complete B-matrix, column by column. 

4.3 Interactive Data 

The prompted inputs are supplied by the user as unformatted, 
list-directed inputs. The following is a list of the prompts 
supplied by the program along v/ith a brief explanation of 
the requested data. 

1) "INPUT TOLERANCE" - The tolerance Is supplied in units 
of the surface definition, e.g., inches or centimeters. 

This number controls the accuracy of each calculated 
point relative to the true intersection curve. 

2) "INPUT SMAX AND TMAX" - SMAX is the maximum allowable 
step size in arc length betv/een the calculated inter- 
section points. This number acts as a constraint on 
the automatic step siiie selector algorithm described 
in Section 3.5. TMAX is the desired slope change 
between consectitive calculated points. In the 
notation of Section 3.5, se = TMAX. 

3) "INPUT NGRID AND NSTEP" - NGPIO establishes the hunting 
grid described in Section 2.1 (in the notation of Section 
2.1, NGRID = Ng). NSTEP establishes the number of steps 
along each curve (in the notation of Section 2.1, 

NSTEP a 

4) "INPUT OUTPUT FLAG (1 GE IGF LE 3)" 

I0F=1 - the output array consists of 

ordered quadruples representing 
the parametric coordinates of the 
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Intersection points. 

I0F“2 - the output array consists of para- 
metric cubic curve fit coefficients 
(neometric form) for the parametric 
coordinates along the intersection • 
curve. The curve segments are 
ordered such that all of the loops 
for surface 1 are given, follov/ed 
by all of the loops for surface 2. 

I0F=3 - the output array consists of ordered 

triples giving the physical coordinates 
along the intersection curve. 


4.4 Output Data 

The output data file consists of ti/o parts. The first card 
image is vfritten using FORMAT (1215) and consists of the 
following integer data 

UW = number of v/ords of output data 

LOOP = number of loops in intersection (_<3) 

NPL(3) = numbt.’ of points in each loop 
flCL(3,2) = number of curve segments for each loop 
for each surface (I0F»2) 

The second portion of the output file consists of the array 
A0(I), lal, MVI. For IOF-1 or 2, the output format is (IH , 
6E12 5). The output array may be decoded as follows 

I0F=1 - Each set of four numbers represents 
a point (u, v, s, t) 

I0F=2 - As mentioned above, the intersection 
curves on surface 1 precede the 
curves on surface 2„ Each segment 
is defined by a group of 8 numbers, 
i.e. uo. V'^, ui, Vi, uj, Vo, uf, vf 
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vmere the primes Indicate derivatives 
and the subscript denotes the value 
of the parametric independent variable. 

I0F=3 - Each set of three numbers represents 
a point in physical space (x, y, z). 


5.0 EXAMPLES 

Two examples were chosen in order to verify the softv/are and to serve 
as benchmark cases. The first example involves a simplified model 
of a wine and a fuse lane. Due to normal aircraft symmetry, only half 
of the configuration was modeled. The first example would then pro- 
duce a single, closed loop intersection curve. In order to test the 
program for multi -loop intersections, the second example was taken 
as a pair of orthogonal circular cylinders of different diameter. 

5,1 Wing-Fuselage Intersection (Example 1) 

The wing-fuselage combination selected for the first example 
is shown in Figure 5.1, The fuselage surface (symmetrical 
half) was modeled as a 4x3 mesh of parametric cubic patches. 

The input file for the fuselage surface is given by Table 5.1. 
The v/ing surface was modeled as a 6x2 mesh of patches and its 
input file is shov'n by Table 5.2, The program v<as run with 
the following set of inputs: TOLERANCE = 0.0001 inches, 

SMAX = 2.0 inches, TMAX = 0.1 radian, MfiRID « 3, NSTEP = 10. 

The output file for I0F=1 is given by Table 5.3. Figure 5.2 
shov/5 the intersection plotted in the parametric plane for 
both surfaces. This example illustrates a rather important 
aspect of the hunting phase, namely which surfaces is designated 
as the primary surface. If the wing v/ere translated vertically 
by a small distance, and if the fuselage were taken as the 
primary with NGRID=3, no intersection points would be found. 

To guard against this possibility and to remove the burden 
of designating one surface as primary, both surfaces are run 
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FIGURE 5.1 - Winq-Body Combination 





through the hunting phase in the prinary role. The surface 
locating the most intersection points is then carried throunh 
as the primary surface for the remainder of the computation. 

5.2 Cylinder - Cylinder Intersection (Example 2) 

This example was selected for its simplicity in illustrating 
the multi -loop capability of the program. Figure 5.3 depicts 
the cylinder geometry. 
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Each cylinder is modeled as an 8x1 mesh of parametric cubic 
patches. The surface input files for CYLIHOERl and CYLIN0ER2 
are given by TABLES 5.4 and 5.5. The program was run using 
the same input parameters as the first example with the output 
file for I0F*1 given in Table 5.6. Figure 5.4 illustrates 
the solution in the parametric plane. 
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